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-^ ' Abstract 

Qi^ ' We consider last-passage percolation models in two dimensions, in which the un- 

(~| . derlying weight distribution has a heavy tail of index a < 2. We prove scaling laws 

^ I and asymptotic distributions, both for the passage times and for the shape of optimal 

paths; these are expressed in terms of a family (indexed by a) of "continuous last-passage 
percolation" models in the unit square. In the extreme case a — (corresponding to a 
distribution with slowly varying tail) the asymptotic distribution of the optimal path can 
^ ' be represented by a random self-similar measure on [0, 1], whose multifractal spectrum 

0^ . we compute. By extending the continuous last-passage percolation model to R'^ we ob- 

tain a heavy-tailed analogue of the Airy process, representing the limit of appropriately 
^b ' scaled vectors of passage times to different points in the plane. We give corresponding 

("^ , results for a directed percolation problem based on a-stable Levy processes, and indicate 

^y ' extensions of the results to higher dimensions. 

^^ 

c^ . 1 Introduction 

a 

Directed last-passage percolation in two dimensions has received much attention in recent 
years. In certain specific cases, for example where the weights at each site are i.i.d. with 
exponential or geometric distribution, very precise scahng laws and asymptotic distributions 
C^ ' are now known, both for the passage times and for the shape of optimal paths (see for example 

|l(illl7l B)). Such cases are closely related to the longest increasing subsequence problem, and 
to Markovian interacting particle systems such as the totally asymmetric exclusion process; 
there are also very close links to random matrix theory, for example to the behaviour of 
the largest eigenvalue of a large matrix drawn from the Gaussian Unitary Ensemble (see for 
example |23] for a survey). 

It is believed that the behaviour proved for the exponential and geometric cases should 
be universal, in that the same scaling laws and asymptotic distributions should occur in the 
last-passage percolation model whose underlying weight distribution is from a much more 
general class (certainly including any distribution with an exponentially decaying tail, and 
maybe also those with sufficiently light polynomial tails). The growth models corresponding 
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to these last-passage percolation problems should belong to the Kardar-Parisi-Zhang (KPZ) 
universality class (see for example [TH|). However, only very limited universality results have 
been proved: for example, conditions under which laws of large numbers for the passage times 
(or "shape theorems" ) hold, and asymptotics for passage-times close to the boundary of the 
quadrant [221 El EI- 

In this paper, we study cases in which the tail of the weight distribution is sufficiently 
heavy that such shape theorems fail, and which certainly fall outside the universality class 
described above. Specifically, we assume that the tail of the weight distribution is regularly 
varying with index a < 2. We describe a family of "continuous last-passage percolation" 
models (indexed by a), and use them to provide scaling laws and asymptotic distributions for 
the discrete models, both for the passage times and for the shape of the optimal paths. Thus 
we have a universality result for these heavy-tailed models as the only information required 
to determine the scaling limits is the parameter a. 

One example of an application where such a heavy-tailed assumption is very natural is in 
the use of last-passage models to represent networks of queues in tandem. The vertex weights 
in the percolation models correspond to service times in the queueing systems, and passage 
times in the percolation models correspond to the total time spent in the queueing system by 
particular customers; see for example |12ll510n|. 

In Section [21 we define the discrete last-passage percolation model precisely; we then 
describe the continuous last-passage model and state our main convergence results. We 
also derive from the continuous model a stationary process which can be seen as a heavy- 
tailed analogue of the Airy process (which was developed by Prahofer and Spohn [2S| and 
Johansson JHDi and which gives a process limit for vectors of passage times to different 
points, appropriately scaled. 

The proofs of the convergence results for passage-times are given in Section [21 and those 
for the optimal paths are given in Sectional The results on the heavy-tailed Airy process are 
proved in section [3 

In Sectioniniwe explore the case where the tail of the weight distribution is slowly varying 
(i.e. Q = 0). It is no longer possible to provide a non-degenerate limiting distribution for 
the passage times; however, asymptotics for the form of the optimal paths are still possible, 
and in fact the distribution of the limiting path that arises can be described in a particularly 
simple and algorithmic way. As the path is increasing it can be thought of as the distribution 
of a random measure on [0,1]; this measure is self-similar and we compute its multifractal 
spectrum. 

The Brownian directed percolation model has recently been much studied in various con- 
texts (see for example [HI ^1 El [^). In Section [5l we discuss a related model in which 
Brownian motion is replaced by an a-stable Levy process, and we again prove distributional 
convergence to the continuous last passage percolation model. 

The bulk of the paper concerns the case of two-dimensional last-passage percolation. 
However, almost all of the results extend easily to d dimensions, d > 3, and now apply 
for a < d. We indicate these extensions in Section [HI 

Simulations of trees of optimal paths, of the limiting path for a — Q, and of heavy-tailed 
Airy processes are given in Sections [21 [HI and [3 



2 Main results 

2.1 Definition of the discrete problem 

Let -F be a distribution function. We will assume that the tail of the distribution F is regularly 
varying with index a S (0, 2); that is, for all i > 0, 

1 - F(tx) 

> t as x ^ oo. 

1 - F{x) 

We will also assume (merely for convenience) that i^ is a continuous distribution and -F(O) — 0. 

The discrete last-passage percolation model with underlying weight distribution F is usu- 
ally defined as follows. 

Let X{i,j),i,j G N be i.i.d. with common distribution F. The quantity X{i,j) represents 
the weight at the site (z, j) G Z^. 

For n G N, we will define the quantity T^"' , the last-passage time between (1,1) and 
{n,n). Let n„ be the set of directed paths between (1, 1) and {n,n). Each such path begins 
at (1, 1) and ends at (n, n), and each step consists of increasing one of the two coordinates by 
1. That is, for any tt G n„, we can write t: — (wi, V2, . . . , f2n), where vi = (1, 1), V2n = {n, n), 
and, for each « = 1, . . . , 2n — 1, u^+i — Vi is either (1, 0) or (0, 1). 

The weight of such a path is the sum of the weights X{i,j) associated with the points 
(i, j) in the path. Then T*^"^ is the maximal weight of a directed path between (1, 1) and 
(ri, n); that is: 

r(") == maxyx(w). (2.1) 

vGtt 

Note that T^") depends only on the weights X(u), u G {1, . . . , n}^. 

2.2 Continuous model 

We start with an alternative representation of the discrete model. Let M^^' > M^ > • • • > 
M 2 be the order statistics, written in decreasing order, from an i.i.d. sample of size n^ from 
the distribution F. 

Consider the set {l/n,2/n, . . . ,{n — l)/ri, 1}^ C [0,1]^, of size n^. Let the sequence 
Y]^^" , . . . , Y 2 consist of a random ordering of the points of this set, chosen uniformly from 
the (n^)! possibilities. We regard F/"^ as the location of the ith largest weight M^' (and we 
have scaled so that all points lie in the box [0, 1]^). 

For two points y,y' G [0,1]^, we say that y and y' are compatible, and write y ~ y', if 
y, y' are partially ordered in that either y < y' co-ordinatewise, or y' < y co-ordinatewisc. 
(Informally, one of y and y' is below and to the left of the other). An increasing path 
will consist of a set of points such that every pair of points in the set is compatible. We 
describe the collection of increasing paths by the collection C'-"\ which depends on the points 
y/"\...,H"^ alone: 

C(n) ^ c(")(n'"\ • • ■ , Y^^^) = {a C {1, . . . , n^} such that for all i,j G A, y/"^ ~ f/"^} . 

(2.2) 



Now we can give a new definition for T'"', equivalent (in distribution) to H2.1|l : 

T(")= max yA//"^ (2.3) 

We formulate the limiting continuous model by defining the distribution of a random 
variable T in an analogous way. 

First let Yi,Y2, . . . be an i.i.d. sequence, with each Yi uniformly distributed on the square 
[0, 1]^. Let Wi, W2, ... be an i.i.d. sequence of exponential random variables with mean 1 
(independent of the (Y,)). Now write, for each fc e N, Mk = (Wi + . . . + WkY^^°'. (Then 
with probability 1, Mfe > Affe+i for each fc, and M^ — * as /c -^- oo). The motivation for this 
definition is given by equation (|2.6I) below. Affe is the fcth largest weight, which we imagine 
positioned at the point Yfe G [0, 1]^. (The set of locations Yjt is of course dense in [0, 1]^ with 
probability 1). 

Analogously to p. 2(1 . we represent the set of increasing paths by the collection C: 

C = C(Yi, 1^2, . . . ) = {A C {1, 2, . . . } such that for all i, j ^ A,Y, ^ Y,} . (2.4) 

Then define 

T=supyAfi. (2.5) 

Remark: Note that one could equivalently define T in (|2.5() as the sup of the weight of finite 
increasing paths A, since either the sup is finite in which case the weight of any infinite 
path can be arbitrarily closely approximated by that of a finite path, or the sup is infinite in 
which case one can find a finite path with an arbitrarily large weight. In particular T can be 
seen as the supremum of a countable family of measurable random variables, and so is itself 
measurable. In Section |3 an equivalent construction of the continuous last-passage problem 
is given using a Poisson random measure approach rather than the sequence of i.i.d. uniform 
positions in the unit square described above. 

2.3 Convergence results 

First note that for all fc, 



V/"),yi"),...ri"))^(yi,r2,...,rfe) 



in distribution, as n ^ 00. 

Now define ajv — _f(^i) {\ — i) . (As an example, if the weight distribution F is Pareto(a), 
with F(x) = 1 — a;~", then ajq = N^^"'. In general, limAr^oo logflAr/logA^ = l/o^)- 

Recall that Af^" is the A:th largest value from a sample of size -n? from the distribution 
F. We write Ml = a~2Ml . Then from classical extreme value theory we have that, for 
all fc, 

(a^i("\ A^^"\ . . . A/f )) -. (Afi, A/2, . . . , Mfc) (2.6) 

in distribution, as n ^ cx) (see for example Section 9.4 of |l()l). 

In particular, M^ is asymptotically of the order of a„2, for any i. (For example, for the 
Pareto distribution F{x) = 1 — x^" mentioned above, we have a„2 — n?''^). Since certainly 



y(") > Mj;"% we have that T^"^ grows asymptoticaUy at least on the order of a„2. In fact, 
we win show that this lower bound gives the right order of magnitude. 



Specifically, let f(") = a;^2^r(") = sup^gc(") E 



i^A 



M, 



(n) 



Then we will show: 



Theorem 2.1 The random variable T defined at \2. 5|) is almost surely finite, and T(") — > T 
in distribution as n ~* oo. 

For comparison, one can consider the case of a lighter tail. If /q°°[1 — F{x)]^^^dx < oo 
(this condition is very slightly stronger than the existence of a finite second moment) then 
a law of large numbers holds: n~^T^'^' — > 7 as n — > oo for some deterministic 7 j21| . If the 
weights are exponential with mean 1, then ^ — A 27 , and then in fact one has the much finer 
convergence result that ri^^/'^(T'^") — An) converges in distribution as n ^ 00, to the GUE 
Tracy- Widom distribution |16| . 
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Figure 2.1: A simulation of the last-passage percolation model with F given by a Pareto 
distribution with index 1. The "tree" consisting of optimal paths from (1, 1) to (i, j), for all 
1 < *ii < 159 is displayed; the thickened path is the optimal path from (1, 1) to (159, 159). 
This represents p(")* in the language of Theorem 14.41 
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Figure 2.2: As Figure I^ITl but now with F given by an exponential distribution. As n grows, 
the distribution of the optimal path from (1, 1) to {n,n) becomes concentrated around the 
diagonal - the deviations of the path from a straight line are on the order of n^^^. This 
contrasts with the deviations on the order of n observed in the heavy-tailed case, where the 
limiting path distribution is non-degenerate (as illustrated in Figure ITT)) . 



We now outline the results on path convergence, which are given in full in Section 0] 
We will show that the optimal path for the continuous model is well defined; that is, with 
probability 1 there exists a unique A* Cz C such that T = X^ieA* ^-^i (attaining the sup in 
(|2.5ll ). Then there is a unique closed connected set P* C [0, 1]^ which contains all the points 
Yi,i G A* and which itself has the directed path property (i.e. y ^ y' for all y,y' € P*). 

Analogously, let A''^^* be the optimal path for the discrete model, achieving the max in 
(|2.3|l (Since the weight distribution F is continuous, the finitely many increasing paths all 
have different weights a.s., and so this maximizing path is a.s. unique). Let p(")* be obtained 
by linear interpolation between the locations {1^ : i € yl(")*} of weights used in this optimal 
path, taken in increasing order. Then we show that the distribution of p(")* converges to 
that of P* as n ^ oo (under the Hausdorff metric on closed subsets of [0, 1]^). 



In the corresponding situation for weights with exponential distribution, the optimal path 
converges instead to a trivial limit, the straight line from (0,0) to (1,1)- In general, the 
deviations of the optimal path from (1, 1) to (n, n) in the discrete model are expected to be of 
the order of n^/^ in cases falling into the KPZ universality class (proved rigorously in certain 
cases ^lE]), rather than on the order of n as we see in the heavy-tailed case. See Figures ITTI 
and 12. 21 for simulations of optimal paths in the cases of weight distributions which are Pareto 
and exponential. 

2.4 Last-passage random fields and the heavy-tailed Airy process 

In Section[71we show how the results described above can be extended to give the multivariate 
convergence of vectors of passage times to different points. 

For x,y > 0, let T^")(x,?/) be the maximal weight of a path from (1,1) to ([nx], [nj/]). 
(So for example the quantity T(") defined in IJOJ is equal to T(")(1, 1)). 

Define also f ("^(x,?/) = a-^T''"\x,y) as before. 

Then we will construct a random field {T{x, y), x,y > 0}, using a Poisson random measure 
construction rather than the sequence of points ordered in decreasing order of weight above, 
in such a way that 

{f(")(x,y),x,y>0}^{r(a;,y),x,2;>0} 

as n —^ cxD, in the sense of convergence of finite-dimensional distributions. 

From a scaling property of the distribution of the weights one has further that the random 
field defined by 

e(w, v) = exp f-^^^\ T (e", e") 
is stationary on M^. The convergence above can be rewritten as 

J exp (-'^^-^\ f(") (e", e") , u, w e M I ^ {e(u, w), u, w € r}. (2.7) 

To remove the multiplicative factor on the LHS of H2.7|l . one can look at a line u+v — const; 
for example, the process Hy — 6(y, —y). This process is stationary in y e M, and we obtain 
the weak convergence 

{a;,ir(") (e^e-^) ,y e m} ^ {e(j/,-j/),y e M}. 

This gives an analogy with the "Airy process" ^Hl [21] i which arises for example in the case 
where the underlying weight distribution is exponential (with mean 1 , say) . There one obtains 
a stationary process limit for the quantities 

1^-1/3 ry(n)(^ _^ yn-l/^ 1 - yn-i/3) - Ai^ , y e »}, 

whose marginals are given by the GUE Tracy- Widoni distribution. Simulations of the "heavy- 
tailed Airy process" Hy are given in Figures [7.117.31 

We also obtain estimates on the moments and correlations of the random field T, showing 
for example that Er(x,y)^ < oo for all /3 < a and giving bounds for E|r(x,y) — T(a;',y')|'^. 



We note that the path convergence described above could also be extended to the mul- 
tivariate setting, to describe the convergence of the distribution of trees of optimal paths to 
a continuous tree structure given by the set of optimal paths in the continuous last-passage 
percolation model. However we do not pursue this further in this paper. 

3 Convergence of the last-passage time distribution 

To establish the convergence in Theorem 12.11 we will work with approximations to T and 
T(") which depend only on the k largest weights. First, define 

Cu = C{Yi,Y2, . . . , Ffe) == {A C {1, 2, . . . , k} such that for all i,j G A, F, - Yj} . 
^(n) ^ c(")(y/"\ . . . , y^(;]^,) = {a C {1, . . . , fc A n"} such that for aU z, j e A, f/") ^ f/")} . 



forC, 



Note that in fact Cfe == {A e C : A C {1, . . . , A:}} = {A n {1, . . . , A:} : A e C}, and similarly 

(n) 
k ■ 

Now let 

Tk = sup Y^ M„ 

^'^'^ i&A,i<k 

and 

ri")= sup Y. Mt^- 

Note that indeed Tk depends only on (Mi, . . . , Mk) and (Fi, . . . , Yfe), while T^"^ depends only 
on (m("\ . . . , m(")) and (f/"), . . . , F,^"'). 

As before, define f^"^ = a-^T^"'' . 

We also define the "remainder terms" Sk, 5^" by 

S-fe^sup Y. ^^« and -^i"^ = s'^P Y. ^^i"^- 
'*^^zeA,z>fc '^^'^^"^ i£A,i>k 

Write also Sl^^ =ajsl^\ 

Lemma 3.1 With probability 1, Sk < oo for all k > 0, and Sk -^ as k -^ oo. 

In particular, putting k — we will have that T < oo a.s. (Later on, we will show more, 
namely that ET^^ < oo for all < /3 < a; see Proposition 17. 2|) . 
We will also have that Tk -^ T a.s. as fc — > cx), since, for all k, 



< T - Tfe = sup ^ Af, - sup Y ^^» 
^^^ieA ^^^ ieA,i<k 

< sup Y ^^i 

Aec ^^ 
= Sk- 



'^'^'^ ieA,i>k 



The convergence in Theorem 12.11 will then follow from the following two results, which 
provide control over Tk — T\ ' and T^"^ — T'"^ for appropriate k: 



Proposition 3.2 Let e > and k be fixed. Then for all n sufficiently large, say n > Nk{e), 
there is a coupling of the continuous model and the discrete model indexed by n under which 



(Elh-^/"1l>^) 



> e < e, 



<e, 



< e. 



Proposition 3.3 Let e > 0. Then for k sufficiently large, 



-An) 



(3.1) 

(3.2) 
(3.3) 



{sr > ^ ) < ^ 



for all n. 



Proof of Lemma 13. II 

First, we define Li = sup^^^ jA n {1, . . . ,i}|. Li is tlie largest number of tlie points 
Yi, . . . ,Yi (the locations of tiie i largest weights) that can be included in an increasing path. 
Note that the collection {Li) is independent of the collection {Mi). Li has the distribution 
of the "longest increasing subsequence" of a random permutation of length i. In particular, 
there is a constant c such that, for all i, KLi < c^fi and EL^ < cz; also, Li/^/i ^ 2 in 
distribution. See for example 1 for a survey. 

We will also write Uk — X^i^fc+i Li{Mi — A/^+i) for each fc > 0. Fix A G C, and define 
Ri = \An{l,.. . , i}\. Then L{i £ A) = R^ - R^-i, and by definition i?j < L^. 

We have 



M,, 



ieA,i>k 



iGA,k<i<n 



i=k+l 



= lim 



y M,= lim y 

n 

= lim V MiL{i e A) 

n—>-oo ^—^^ 
n 

= lim V Ah{R,-R,^i) 
=fe+i 

n-l 

-Mk+lRk+ y Rr{M, - M,+i) + MnRn 

i=k+l 
•i-l 

< lim y Ri{M,-M,+i)+liminiMnRn 

i=fc+l 
n-l 

< lim y Li(Mi -Mi+i) + liminfM„L„ 

i=k+l 

= Uk + liminf M„L„. 



Now lim infn^oo Afn^n = a.s.; this follows, for example, since (by the law of large 
numbers) M„ ~ n^^'" a.s. (with a < 2), and since Ln/y/n converges in distribution to a 
constant. 



Since the inequality above holds for any A E C, we therefore have that Sk < Uk a.s., for 
any k. To conclude the proof, we will show that with probability 1, Uk < oo for all k, and 
C/fc — > as fc ^ oo. (In fact, as soon as f/i < oo, we necessarily have that C/fc — > as fc — > oo, 
since the quantity Uk is the "remainder" from index fc + 1 onwards in the infinite sum Ui ; if 
the infinite sum Ui converges, then by definition these remainders tend to 0). 

Hence it's enough that Uk < oo a.s., for all k. Specifically, we'll show that KUk is finite 
whenever k > 1/a. Then certainly Uk < oo a.s. for such fc, and in fact Ur < oo a.s. for all r, 
since if r < /c, Ur — Uk is the sum of only finitely many terms. 

By independence of the collections {Li) and (Mi), 



EUk= J2 EL,(EM,-EM,+i) 

OO 

< Y^ ci^^^{EMi-EM,+i). (3.4) 

i=k+l 

Now Mr has the distribution of {Vr)~^^°', where Vr has Gamma(r, 1) distribution. We 
then obtain 



/.oo -, 

EM,. = / _-x;'^-ie-%-i/"dw 

Jo r r) 



/o r(r) 

= r(r-l/a)/r(r). (3.5) 

Using the identity T{z + 1) = zT{z) and the fact that the gamma function is log convex, one 
has 

(^ - 1)" < ^^^^7^ < (x + a)" (3.6) 

for a; > 1, a < 0. Then 

E M,. - E Mr+i = 



T{x) - 


yu, -r u,/ 




Tir-l/a) 


r(r + l-l/a) 


r(r) 

r(r-l/a) 


r(r4 

r 


-1) 


r( 
1 r(r-i/ 


r) 
a) 





ar r(r) 



< 

ar 



-(r-l/a-l)-!/". 



Returning to 1)3. 4|). we have 



i=k+l 



which is finite for all k > 1/a (since a < 2). D 

Proof of Theorem[2IT] Wc will find a couphng of f (") and T for each n such that f^'^)~T -> 
in probability as n —^ oo. 

For each n S N, define A:„ — max{/c : n > Nk{l/k)}. 
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Then fc„ — > oo as n — > oo, and, for all n, n > Nk„{l/kn). Hence from Propositions 13.21 
and 13.31 and from Lemma l3. II there are couplings such that, as n ^ oo, 



Eh 

i=l 

fen 

EN 



M. 



(") 







(3.7) 



K 



(") 



5; 



(n) 



k„ 
'5'fe„ 



in probability, and 



Now 



(") 



T - T(") = (T - TfeJ + (Tfe„ - T^^ 

( 



0. 



j,(ri) _ rp{n) 



We have |T - Tfc„| < Sk„ and |f^"^ - f ("^| < 5*^"^ so to show that the LHS converges to 



in probability as desired, it remains to show that I T^^^ — T^ 
We have 

Tk^ = max y Mi and 



(«) 



max 



rf{n) 

fe„ 



ieA 



max 

(") 



Aec) 



-> in probability. 

E^i 



fin) 



SO if Ct„ 



Ci"\ then 



Tk^ - t: 



(") 



Since P Cfc„ ^ C 



"fen 



and Y.rli 



M,; - M 



(«) I . 



(n) 



M. 



in probability, we are done. 



D 



To complete the proof it remains to prove Propositions 13.21 and 13.31 

3.1 Convergence of f;^"^ to T^ 

Proof of Proposition IX^ 

We have (m}"\ . . . , M^"\ r/"\ . . . , r^l"^) -^ (Afi, . . . , Mfc, Yi, . . . , n) in distribution as 
n -^ 00. By the Skorohod Representation Theorem, we can define all the variables on the 
same space in such a way that the convergence occurs almost surely. Then indeed (|3.1|) and 
(|3.2|l must hold for large enough n. 

Note that since the variables Yi are i.i.d. uniform on [0, 1]^, there are almost surely no two 
i and j such that Yi{d) = Yj{d) for d = 1 or 2. 

Thus if we perturb the point (Yi, . . . , Yfe) by a small enough amount, the orderings of all 
the coordinates remain the same, and the set Ck of increasing paths is unchanged. In fact, if 

Z l<i,j<k,i^j a— 1,2 

-f Yi a.s. on the joint probability space for all 1 < i < k, 
we then have C^"^ = Ck eventually, with probability 1. Thus (|3.3|l must also hold for all large 
enough n, as desired. D 



max ||r,-y/")| 

l<i<fe * 

then Cj! — Ck- Since we have F/" - 
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3.2 Convergence of T^") - T^^"^ to 

Our aim in this section is to prove Proposition |^2I 
Define the "good event" S^" : 

Sl"^ = I F-i d-^] < Afi") < F-i (l - X] for aU fc < r < nn . (3.8) 



Lemma 3.4 P [Bj!^ ) ^ ^ as k ^ oo, uniformly in n. 
Proof: 

P (ei") fails ) < P (Mi% > F-i 1 - 1 ) ) + ^ P M(") < ^-M 1 - I 

L"'/2J 

= P (Binomial(n^ l/n^) > k + l) + ^ P (Binomial(n^ 2r/n2) < r) 

r=fc+l 

using Markov's inequahty for the first term and an estimate from Corollary 2.3 of ^H] for the 
second. The RHS tends to as fc — > oo, uniformly in n, as required. D 

Now, we will prove a bound on the expectation of Sj. in terms of the "order statistics" 
Mr^"\ Define 

l|"^ = max \An{l,2,...,i}\. 

Recall that F^^"^ e {1, 2,...,n}^ is the location of the rth largest weight, Mi"\ Thus, L,^"'' 
is the maximum number of the points Y^^^ , ■ ■ ■ , Yi that can be included in an increasing 
path. 

Note that the collection (Lr ■^)i<r<„2 is a function of the values ly" alone; in particular 

(n)^ (n) 

it is independent of the weights M,- and of the events Br . 

Lemma 3.5 There is a constant c independent ofm andn such thatELm < c^/m, whenever 
1 < m < n^ . 

Proof: The distribution of {^^ , . . . , Ym } is uniform over the subsets of {1, ... , n}^ of size 
m, and Lr is the maximum number of the points i^" , . . . , ly" that can be included in an 
increasing path. 

We compare this with the last-passage percolation problem in {1, . . . ,n}'^ with i.i.d. 
Bernoulli(p) weights. 

We have the following representation for the expectation of the passage time for such a 
problem: 

E Bcr(p)T(n, n)^Y. ^(^p,"^ = ^)^ ^^"^ ' (3-9) 

r 

where -Bp.„2 has Binomial(p, n^) distribution, since, conditional on the event that exactly r 
of the n^ weights have value 1, the set of positions of the weights with value 1 is uniformly 
distributed among all the subsets of {1, ... , n}^ of size r. 
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Proposition 2.2 of TT shows that there exists a constant ci such that 

IEBcr(p)T(n,n)<cipi/2^ (3.10) 

for aU n, p. 

Given 1 < m < n^, set p = mm{2m/n'^, 1). Using l|0)l . (|XTU|) and the fact that Eil"^ is 
increasing in r, we obtain 

P(i?p,„. > m)EL(;^) < ^P(Bp,„2 = r)E4") 

r 

1 /? 

< ci (niin(2?n/n , 1)) n 

< C2\fm. 

To complete the proof, it then sufhces to bound P(_Bp „2 > rn) away from uniformly in 

1 < 77T. < n^. 

If m > n2/2 then p = 1 and P(Sp,„2 > to) = 1. 

For 1 < TO, < r? jl., we have pv? 12 = m, and we use the estimate 

P(Bp,„2 < pnV2) < exp(-pnV8), 

(see for example Theorem 2.1 of QH), to give P(i?p^„2 > to,) > 1 — exp(— 1/8) uniformly in 
1 < TO, < n^/2 as desired. D 

Lemma 3.6 



r=k+2 

Proof: The argument is similar to the proof of Lemma l3. II 

Let A achieve the max in the definition of Sj^ , so that 5^^^" — J2i£A i>k ^i ■ 

Define i?,; — |A n {1, 2, . . . , i}\ for each i. Then Ri — Ri-i = I(i G A), and by definition 

Ri < Li for each i. We then have 



5f'= Y. ^^i 

iGA,i>k 
i=k+l 









i=fc+l 

Since i?|") < l|") and M^ ^ > M^l\. 

We now take expectations, restricted to the event B)^\ using Lemma [3.51 and the inde- 
pendence of the Lr from the Mr '■ 



{si"^;B^:^) < Y: IEL|"^ [e (M(");6i")) -E (m^;!;^!' 



E 



n) 
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i=k+l 



E (Mf"^;4") 



E(M(;1;Si"^ 



= c^/fcTTE(Mi:;\;Si") 



n 



E 



i=k+2 



lEii'J^E (Mf^"^;4"^ 



r("). R(n) 



n 



<cvmE (A?(:;\;sfn +c ^ *-i/2ie (Mr';S, 



«■>(«). J^(n) 



i=*:+2 



since i^^^ ^ {i ~ 1)^^^ < i ^^^. This is the required result. 



D 



The next lemma gives an estimate on the tail behaviour of the weight distribution using 
the regular variation condition. Note that if the weight distribution were Pareto(Q;), then 
F-^{u) = (1 - u)"!/", and one then has exactly F-^{ui) = F-^{ua) [(1 - ui)/(l - uq)]"^/" 
for any uq, ui. 

Lemma 3.7 For any S > 0, there exists U{6) < 1 such that for all Uq, Ui with U{S) < Ui < 
uo, 



F-\ui) < 2F-\uo) 



1 — Ui 



'^+S 



, 1 - "0 . 

Proof: Fix s > 1 sufficiently small that s^/"^"^ < 2. 

From the fact that the tail of F is regularly varying with index a, the following property 
holds: if ui is sufficiently close to 1 (at least U{S), say), then for all (1 — uq) < (1 — ui)/s, 

Iterating, one obtains that if U{S) < ui < uq, then 



F-^{uo) 



<^s- 
= exp 
< exp 



Llog,((l-«i)/(l-«o))J 



+S] (logs) 



log((l~«i)/(l-«o)) 



logs 



< 2 



1 — Ml 

1 -Uo 

1 — Ul 

1 -Uo 



1 \ / 1 - Ui 

\-S\ log- logs 

a J \ 1 - Uo 



S' 



as required. 



D 



!>(») 



Finally we use this tail estimate to control the expectation of the variables Mr , restricted 
to the "good set" b'j^K 

Lemma 3.8 Let 6 > 0. Then there exist cq, ci and C2 > such that 
E (Mi");^^"^) < cor-i+' + c^ajlir > c^r?) 
for all n, k, r satisfying 2(1 + 1/a) < k < r <n^ . 
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Proof: Let U{S) be as in Lemnia l3.7l and set C2 == {l-U{S))/2. Then r < 02^^ <=> l-2i^/n'^ > 
U{S). 

Note that if 

max i C/((5), l-^i<u<l-^, 

then, by Lemma |3. 71 

< 2a„2n-2/"(l - M)-i/"(2r)*. (3.11) 

Since r > k, we have that 

«i'"c{f-'(i-|)<M,<".<f-(i--i; 

Hence 



e(m^");6^") 



''-^fc 



< a-J^ ( Af(");F-i ( 1 - J) < M^-^ < f "^ (l - ^ ) ) (3.12) 



.^e(f- ([/(")) ;l- 5 <C/(")<1--L 



-1 Z"'"'^"' -1 

= a„2 / i^ (u)/r.„2(u)dw 

= a-Jl 1 - — < U{5) F-\U{5)) + a~l / i^-i(w)/,;„2(«)du 

I " J Jmax{L'(<5)a-2r/n2} 

.l-l/n" 

Jmax{U{S)S-2r/n^} 

< cil {r > cin^) + C3a„2n-2/"/ / (1 - M)"^/"/^;„2(M)du, (3.13) 

Jo 

where ci = F~^([/(5)) and C3 = 2^+'^, and where jr;n? is the density function of the rth 
largest from an i.i.d. sample of size r? from the uniform distribution on [0, 1]. 
Now 



and, since r — 1 — 1/a > 0, one then has 

Jo r(n^ - r + l)r(r) Jq 



r(n2 + 1) ^ r(n2 + 1 - 1/a) 

(3.14) 



r(n2 - r + l)r(r) Vr(7i2 - r + l)r(r - 1/a) 
r(n2 + i) r(r-l/a) 



r(n2 + 1 - 1/a) r(r) 
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Using (|?^ . the RHS of (|m|) is bounded above by 

[n^ + 1)1/" 
(r- l-l/a)i/"' 

Since r > 2(1 + 1/a), this is in turn no greater than (4n^/r) . Inserting this into H3.13(l 
gives the desired result. D 

Proof of Proposition 1^31 We may assume that k <n?, since if fc > n^ then S^ — 0. 
Fix e > 0, and fix some 5 < ■^~ \- Using Markov's inequality, we have 

P (^I"^ > e) < F(6(") fails ) + e-^E (^^"^Si"^) . 

By Lemma 13.41 the first term tends to as fc ^ oo, uniformly in n. For the second term, 
Lemmas I3.m and l3.8l combine to give 



r=k+2 

+ cci{k + 1)^ a:^2 + ccia:^2 ^ r"^. 

From the choice of 5 and the fact that a„2/n ^ oo as n — > cx), one obtains that all four terms 
on the RHS tend to as fc — > oo uniformly in n such that k <ii? . 

Hence indeed P ( S^^ > e ) < e for all large enough k, uniformly in n such that fc < n^, as 
required. D 



4 Path convergence 

In this section we will state and prove results describing the convergence of the distribution 
of the optimal paths for the discrete models to that for the limiting continuous model. 
First we note that the optimal path for the continuous model is well-defined: 

Proposition 4.1 With probability 1, there exists a unique A* ^ C such that T = X^ieA* -^*- 
We will also define yl^")* as the set that achieves the maximum in 



T(")= max^^Aff"^ 



Aec(") "^' 



(or in the equivalent expression for T*-"'). (Since the weight distribution is assumed to be 
continuous, this optimal set is almost surely unique). 

It will be useful to extend the sequences {Y^ )i and {M^ )i to all i G N (rather than 
only i < n^); for example, we can put Mj = 0, Y^ = (0, 0) for all i >n'^. We will consider 
always the product topology when looking at convergence of such infinite sequences. 

We also use the product topology on S, the set of subsets of N. Thus, given A and a 
sequence A^ in S, we have A^ -^ A if, for every m, A^ H {1, . . . , ni} is equal to An {1, ... , m} 
for all large enough fc. One has easily that any sequence A^ has at least one limit point, 
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and also that if Ak G C for each k then every limit point is also in C (since if the limit point 
contains i and j, then i and j are in Ak for some fc, and hence Yi ^ Yj). 

The following theorem is our first path convergence result. Later (in Theorem I4.4|l we 
will use it to prove a more direct convergence result concerning the optimal paths viewed as 
random subsets of [0, 1]^. 

Theorem 4.2 ('(r/"^)^gpj, A^")*") -> ({Yi)^^j^^,AA in distribution as n ^ oo. 
Before proving Proposition 14. II and Theorem 14. 21 we need the following fact: 

Lemma 4.3 With probability 1, the following holds: if Aj is a sequence in C converging to a 
limit A, then limj^oo Y^ieAj ^'^i = Z^ieA ^i- 

[N.B. this result is not true in general for sequences Aj in S (unless a < 1 so that Yl ^'h < '^ 

a.s.)] 

Proof: If Aj n {1, . . . , m} = A n {1, . . . , to} then 



i£Ai iGA 



< sup yj Mi = S„ 



^'^'^ ieA,i>m 



But S„i — > a.s. as m -^ oo (from Lemma I3.1|l . and Aj C\ {1, . . . , to,} = A n {1, . . . , ttt.} 
eventually for all to, so we are done. D 

Proof of Proposition |4. II Recall that 

Tk == sup ^ Mj. 

Since the sum on the RHS depends only on the intersection of A with {1, . . . , /c}, we only 
need to consider the max over finitely many A C {1, . . . , fc}. Thus there exists some A^ which 
achieves the sup. 

We consider the sequence A^.. As observed above, this sequence has at least one limit 
point A* e C, and by LemmaOl HieA* ^^» = limfe^oo Y^ieA' ^'^i = l™fc^oo Tk = T. 

Now we wish to show that in fact a unique A* achieves the sum T. 

Suppose instead that there are two such optimising sets in C. Then there is some k that 
is contained in one but not the other. 

In that case, 

sup 2_, ^^i — sup yj -^^i 



= Mk+ sup J2 ^*' 



Aec,keA.^^^_^i_, 



which gives 



Mk= sup VMi- sup V M,. (4.1) 

We will show that this event has probability for each fc; then, by countable additivity, we 
are done. 
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The RHS of H4.1|l does not depend on M^; in fact, it is a function of the collection 
((Mi)jgj>[y{j,j, (Ki)igN)- We condition on the value of this collection. Then the RHS is a 
constant, while the random variable Mk on the LHS has a continuous distribution. (Specif- 
ically, the distribution of M^" conditional on this collection is uniform on the interval 
i^J^k-i^^'Ik+i)^ since the sequence (M^^") forms the points of a Poisson process). Hence 
the event (14.1(1 has probability 0, as required. D 

Proof of Theorem 14.21 We will use the same couplings as in the proof of Theorem 12. II 
As at (13.711 ■ we then have that P (C[" ^ Ck„ ) ^ as n ^ cx3, and that all of the quantities 



i=l 



M, - M, 



(n) 



1=1 



Y,,~Y, 



(n) 



, Sf^ and Sk^ converge to in probability a.s n -^ oo. 



To prove Theorem 14.21 it will then suffice to show in addition that for any to, 
P (a(")* n {1, . . . , m} ^ a* n {1, . . . , ?7i}) ^ 
as n -^- c». For this, it's in turn enough to show that for all r, 

P(re A*,r^ A(")*) ^0 (4.2) 

P(r^ A*,rG A(")*) ->0 (4.3) 

as n -^ oo. We will show (|4.2|) : an analogous argument gives (I4.3|l . 

Define T^^r) = ^^VAeCr^A E,eA ^^»- 

Suppose r e A*. Then r(_r) < T strictly. (Otherwise, there is a sequence of members of 
C, none of which contain r, whose weight converges to T; then (by Lemma IO|) this sequence 
has some limit point, itself a member of C not containing r, which attains the weight T. But 
this contradicts the uniqueness of A* established in ProDOsition l4.1(l . 

IfC^"^=Cfc„, then 



max y^ M 

^ teA 



(") 



< 



max 

Aeci'^\r(^A 



E^I. 



>(") 



si:^ 



ieA 



< max y^ M, 



i>(") 



i£A 



<^(-)+E 



i=l 



[Mt^ 



E 

1=1 



M, 



[m^") 



- Mr 



r.(«) 



s^ 



(4.4) 



and similarly 



max 

A<£C'.">,r(£A 



i£A 



(n) 



> T 



k„ 

E 

z=l 



M,, - M 



(«) 



-Sk„- 



If also r^ A(")*, then 



max } M > max 

^ i6A leA 



E^"'' 



and using (|4.4|l and (|4.5|l we get 

T 



k„ 



Ti-r) <Er^^ 



M, 



(") 



s, 



k„ 



S 






(4.5) 
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So altogether we obtain 

< p (cf^ ^ Cfe) + p < T - r(_,) < ^ |m. 



M^"^ 



^fe. + ^i!^ 



i=l 



We have already observed above that the first probability on the RHS tends to as n — > oo. 
The same is true for the second probability on the RHS, since all the terms on the right of 
the inequality converge to in probability as n ^ oo, while the term in the middle of the 
inequality docs not depend on n. Hence F (r £ A*,r ^ ^(")*^ ^ as n ^ oo, as required. D 

We now turn to the convergence of the paths regarded as subsets of [0, 1]^. 

First let U* = {JieA- ^* ^ {(0, 0), (1, 1)}, and take its closure U*. 

We expect that U* is connected with probability 1, but we don't have a proof. To work 
around this, wc will use the fact that, at least, there is a.s. a unique way to extend U* to 
a connected set while preserving the increasing path property. (Here the increasing path 
property of a set means that if y and y' are two elements of the set then y ^ y'). 

To see this, first note that if U* does "contain jumps" , then none of these jumps can span 
a rectangle of non-zero area. That is, with probability 1 there is no rectangle R of non-zero 
area such that y ^ y' for all y E U* and y' E R. 

For if there were, then R would certainly contain some points Yj] such j could be added 
to A* , increasing the weight of the path by Mj ; this contradicts the maximality of A* . 

So any jumps in U* consist only of horizontal or vertical line segments. These segments 
can all be added to U* while still preserving the increasing path property, and this gives a 
connected set. Conversely, any connected increasing set containing U* must "fill in" these 
jumps. 

Thus, define P* by setting y E P* if: 

(i) y E U\ or 
(ii) there exists y' ,y" E U* with either 

(a) y'{l) = y{l) = y"(l), y'(2) < y(2) < y"(2), or 

(b) y'{2)^yi2)^y"i2),y'il)<y{l)<y"{l). 

This set P* (which wc conjecture to be equal to the closure of IJjg^* Yi w. p. 1) provides 
the distributional limit wc need. 

For each n, we define the object representing the optimal path in the discrete problem 
indexed by n as follows: order the points {Y^ ,i E A'"'*} in increasing order and join 
successive points by a straight line (horizontal or vertical, of length 1/n). Call the resulting 
path P(")*. 

P* and p(")* are regarded as subsets of [0, 1]^ and we use the Hausdorff metric: 

dH{Pi,P2) ^ sup inf \x-y\+ sup inf \x - y\. 
xePi y^P^ xeP2 y^Pi 

Theorem 4.4 p(")* -^ P* in distribution as n ^ oo. 
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Proof: Choose a probability space on which the convergence in Theorem 14.21 occurs almost 
surely. We will show that in this case p(")* — > P* a.s. also. 

First consider any point y G P* . We will show that for all sufficiently large n there is a 
point of p("')* within distance e/2 of y. 

There are two cases to consider. 

First, suppose that y is a limit of some sequence of points Yi G A*. Choose some Yi which 
is with distance e/4. For large enough n, we have 1^/" G A^"'* and \Y^ —Yi\< e/4 and we 
are done. 

Otherwise, y is on a vertical or horizontal line between two points that are limits of 
sequences Y^ €z A*. Call these endpoints y~ and y^. For large enough n, there are i~ and i+ 
such that Yi- , Yi+ e P(")* with |r/"^ -y^\ < e/2 and |r/+"^ - y+| < e/2, as above. Then by 
the increasing path property, the subpath of p(")* joining Y^- and Yi+ passes within e/2 of 
every point on the line segment joining y^ to y+, and hence in particular within e/2 of y as 
required. 

Now for some rn G N set e = 1/m and consider an increasing sequence of points (0, 0) == 
y(o)^y(i)^ . . . jy^^™) — (1^ 1) £ p* such that the Li distance dl(y^-'^ y^-'"*'^^) between successive 
points is exactly e for all j. (This is possible since P* is an increasing path and connected). 

Now for large enough n there is a sequence y^°\y^-^\ . . . ,y(2™) of points of p(")* such 
that di{y^^\y^^'') < e/2 for all j. Then necessarily y("',y(^\ . . . ,y(2™) is itself an increasing 
sequence and rfi(y'^-'\y*^-'+^^) < 2e for all j. 

Using the increasing path property for P*, we have that every point of P* is within Pi 
distance e/2 of one of the y'--^K Then since each y^^^ is within Pi distance e/2 of a point of 
p{n)*^ we have that every point of P* is within e of p(")*. 

Similarly, the increasing path property for p(")* gives that every point of p(")* is within 
ii distance e of one of the y'-'\ and thus in turn within distance 3e/2 of some point of P*. 

Hence d/f (P*, P*^"^*) < 5e/2, for all large enough n. This works for any e = 1/m, so 
p(")* _^ p* as required. D 



5 Stable process directed percolation 

In this section we consider a directed last-passage percolation model based on stable Levy 
processes. This is the stable version of the Brownian directed percolation problem considered 
in piin. 

For n e N, t > 0, consider the random variable 



Lin,t)= sup X!'^*i-i 



0=ta<ti<...<tn=t _ 

where 5*^ = SI — SI, and S^ are i.i.d. a-stable processes for some a £ (0,2). The Brownian 
version of this problem, in which the stable processes are replaced by Brownian motions, gives 
a representation for the largest eigenvalue process in "Hermitian Brownian motion" (a matrix- 
valued process whose marginal at any fixed time has the CUE distribution) and has been much 
studied in various contexts (see for example [5l ll3[[nil24p . We do not have a random matrix 
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interpretation of this stable process version; however, an interesting connection could be to 
the case of Wigner random matrices with heavy-tailed entries considered by Soshnikov in [2Sj , 
where a scaling is obtained for the largest eigenvalues which corresponds to the one we have 
observed for the heavy-tailed last-passage percolation problem. 

We will show that the asymptotic behaviour of the distribution of L{n,t), as n becomes 
large, is again described by our continuous heavy-tailed last-passage directed percolation 
problem. Note that by scaling L{n,tn) = t^'"L{n^n) in distribution and hence we can just 
consider L(n, n). 

The processes S^ have jump measure c+a;~"~^/:r>o +c_|a;|^"^^/:r<o, for some c+ > Q and 
c_ > 0. The jumps play the role of weights for the percolation problem. 

Theorem 5.1 

1/ 



( — J n-'^''^L{n,n)^T 



in distribution as n —> oo, where T is the last-passage time in the continuous last-passage 
percolation model with index a, defined at 1^2.5]) . 

A short argument is available in the case a < 1 (making use of the fact that the sum of 
all positive weights is finite), and we give this first. 

5.1 Case a < 1 

Let M^\ M2 , ... be the set of positive jumps of the processes 5'^, 5^, . . . , S"" on the interval 
[0,n], written in descending order. 

From the form of the jump measure, we can regard the ordered sequence of jumps as a 
Poisson random measure. Thus by a suitable transformation we can write the sequence of 
jumps in terms of a Poisson process and have that for any n, 



l/c 



n-2/"(M(")^Mf\...,Mf\...) 

^ (VF-^/", (PFi + W^r^'^. . . . , (W^i + • • • + W^fc)-!/", . . .), 

where Wi are i.i.d. exponential random variables with mean 1. 

Now let L^" be the maximal weight of a path, if one ignores all the weights except the 
k largest positive weights. Just as in the discrete case, one can show that 



l/a 

n" 



^ ^ -2/a 7-(")+ , rri 



c • ^ 

in distribution as n — > 00, and one also has T\^ —* T \vl distribution as /c ^ c», where Tj, and 
T are the last passage times for the continuous problem as defined before. 

Now let L(")+ be the maximal weight of a path, if one considers all the positive weights 
but ignores all the negative ones. Then 

00 

■r=fe+l 
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Now the distribution of n ^' " X^^fc+i Mr does not depend on n, and converges to in 
distribution as fc ^ oo (since the sum of all the positive weights is a.s. finite). So we have 



1/q 



in distribution. 

Now consider the optimal path attaining _L'^"^+. Consider the sum of (the absolute values 
of) all the negative weights along the path; call it S"*^"^". Since the positive and negative 
weights occur independently, S"^"^" has just the same distribution as the sum of the negative 
weights for a single stable process between times and n. This is finite and on the scale n^'" 
(in fact, the distribution of n^^/^S"*^"'^ is independent of n). So certainly n^^/^S"*^"'^ -^ 
in distribution as n ^ oo. Since L(")+ — S'*^"'^ < L{n,n) < L(")+, we obtain 

1/ 



( — J n-^^''L{n,n)->T 



in distribution as n ^ oo, as required. 



5.2 Case l<a<2 



Lower bound: 
As before, 



in distribution as 



l/a 

n 



-2/.^(n)+ ^ y^ 



Now consider a path realising L^ in this way; (for definitcness, say the first such path 
in the lexicographic order). 

Let T^" be the total weight of this path, including all weights, and let 5^ = T^ — L)^ . 

The distribution of 5* can be described as follows. Generate the n independent processes 
S^ , . . . , S"' from time to time n. Remove the k largest positive jumps that occur in the n 
processes in time [0,n]. Then S has the distribution of the altered value of S}^, after the k 
largest jumps from the set of processes have been removed. 

However, as n — > oo, the probability that any of the k largest jumps occur in the process 
Si tends to (it is no larger than k/n); so with high probability, the procedure in the 
previous paragraph does not alter the value of 5';'^. Thus the limit in distribution of n^^'^S 
is the distribution of n~^'°'S^ (which is independent of n). In particular, n^^^°'S ^ in 
probability. Thus, for any fc, 

1/q / \ l/a / \ l/a 



in distribution, as n ^ oo. Since L{n, n) > Tj^"" and Tk ^ T a.s k ^ oo, this establishes that 
T is a lower bound for the limit in distribution of I -^ ) n^^'"L(n, n). 
Upper bound: 
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We first need a lemnia on the tail behaviour of the difference between the supremum and 
infimuni of the stable process on an interval: 

Lemma 5.2 Let S be a stable process with index a and jump measure c+x'^"^^Ix>o + 
c^\x\~"~^Ix<o, where c_(- > and c_ > 0. Then 

Pf sup {St - Ss} > x) - ^x"" 

as X ^ oo. That is, the quantity supQ<^,<j<]^{S'( — Ss} has the same positive tail behaviour as 
the size of the largest positive jump of S in [0, 1] (which is also the same as the upper tail of 
Si and o/supo<t<iS'J. 

Proof: Let the running infimum and supremum processes be denoted by It — infjS's : < 
s <t} and Sf = snp{Ss ■ < s < t} respectively. Consider the reflected process Xt = St — It- 
Our aim is to determine the tail behaviour of X^ — sup{Xt : < t < 1}. By Bertoin 
[7] Prop VI. 3, we know that for fixed t, the distribution of Xt is the same as that of Sj^ . By 
standard results (e.g. 7^ Prop VIII.4) we have 

P{Xi >x)^ F{S*i >x)^ ^x-", (5.1) 

a 

as X — > oo. 

Thus we just need to show that X'l has the same tail as Xi. The proof is analogous to 
that of 13 Prop. VIII.4, and we reproduce the argument here. An easy consequence of (|5.1f) 
is that 

liminfP(Xi* >x)x°' > —. 

X — *oo Oi 

Now fix e > and note that the reflected process is Markov by |7j Prop VI. 1. As the stable 
process scales in that Sxt — X^'^St, this property will be inherited by the infimum and hence 
the reflected process itself, giving X\t = \^^°'Xt. Since St — Xt = It is decreasing in t, we 
also have that Xi — Xt > Si — St for all t < 1. Applying these properties and denoting by t^ 
the first hitting time of the interval {x, oo) by the reflected process X, we have 

¥{Xi > (1 - e)x) > P(X* > x,Xi > (1 - e)x) 

> / P(r^ e dt)¥{Xi -Xt> -ex) 
Jo 

> / P{tx e dt)F{Si -St> -ex) 
Jo 

(-1 

F{tx e dt)F{Si-t > -ex) 

10 

> / F{tx e dt)F{Si > -ex) 
Jo 

= F{Xl > x)F{Si > -ex). 
As P(5i > —ex) ^ 1 as X ^ oo we have that 

limsupP(X* > x)x°' < (1 - e)-"c+/a, 



23 



and, as e is arbitrary, we have the result. D 

Now for 1 < i, j < n, define 

X(^,J)= sup {Sl~Sl}. 

j~l<s<t<3 

Let T{n,n) be the maximal passage time for the discrete model with weights X{i,j). Then 
one can see by a direct sample path comparison that i(n, n) < T{n, n). 

Applying Lemma O a^ ^ (-^) N'^'", where a^ = inf{x : T{X{i,j) > x < l/N}. 
Thus, applying Theorem 12. II for the discrete model. 



(^)'^"n-2/"T(n,n)^r 



in distribution as n — > oo. This gives the required upper bound in distribution for the limit 
of L{n, n). 



6 The case a = 0: convergence to the greedy path 

In this section we consider the discrete last-passage percolation model in the case a = 0. The 
distribution F is said to have a slowly varying tail: for all t > 0, 



equivalently, for all s < 1, 



1 - F{tx) 
1 - F{x) ~ 

F-^{l-sv) 



I as X ^ oo; 



oo as w J, 0. (6.1) 



F-^{l~v] 

Now it is no longer possible to find a non-degenerate limit in distribution for T^"' as we 
did in Theorem l2.1l Let M{" be the maximum of an i.i.d. sample from F of size ri^ as before. 
Let bn be any sequence of constants. Then any limit point in distribution of the sequence 
b~^M{' must be concentrated on the set {0, oo}. From Proposition 16 . 21 below, the same is 
true if we replace m["'^ by T^"). 

However, convergence in distribution of the optimal paths can still be obtained. In fact, 
the form of the limiting distribution has a particularly simple description in terms of a "greedy 
algorithm" . We give a multifractal analysis of this limiting object in Section 16.21 

The limiting object is defined as follows. Given the locations Yi,l2, • • • i.i.d. uniform on 
[0,1]^, let C, the set of increasing paths, be defined as at H2.4|) . We now define the greedy 
path A* — A*{Yi,Y2, ■ ■ ■) G C recursively as follows. Let 1 e A* always, and then, given 
A* n {1, . . . ,r}, let r + 1 e A* if and only if Yr+i ^ Yi for every i e A*, i < r. One can 
describe A* as the first member of C in the lexicographic order. 

The discrete problem is defined in terms of the locations iY^') and weights {M^') as 
before. Write Af^"^* for the optimal path as in Section 0] The following theorem gives the 
convergence of these optimal paths to the greedy path: 



Theorem 6.1 



((l^/"^)^^pj,A(")*) -> ((y,)^^p,,^*) m distribution as n 
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Exactly as in Section IT!^ one can also define P* and p(")* to represent the optimal paths 
regarded as subsets of [0,1]^, and obtain the convergence in distribution of p(")* to P* as 
n —^ oo (under the Hausdorff metric). In fact, the situation is considerably simpler here; one 
can simply define P* to be the closure of Ujga* ^*' since by the results of Section |^1 this set 
is connected w.p. 1. 

Theorem 16. II will follow from the next proposition: 

Proposition 6.2 For all r, 



i—r+l 



Proof of Theorem 16.11 Fix some e > 0. As in Proposition 13. 21 for n large enough, we can 
find a coupling such that, with probability at least 1 — e, 



< e 



II ' 

and 

Cf ^ = Ck (6.2) 

Suppose that ^^ holds and also that, for aU r < fc, M,^"'' > E"=r+i ^i"'- Then indeed 

^(")* n{l,...,k}^A*n{l,...,k}, (6.3) 



where A* is the "greedy path". Then using Proposition 16.21 we can find n such that 
holds with probability at least 1 — 2e. This gives the convergence in distribution in Theorem 
O □ 

6.1 Proof of Proposition 16.21 

Lemma 6.3 Let e > and C > 0. There exists U{C, e) < 1 such that if U{C, e) < ui < 
(1 + e)u2 — e, then 

F-l(ui) < F-\u2)il - Ui)-^\l - U2f. 

Proof: Let i = (1 + ef^'. From H6.1|l there exists V > such that for aU v <V, 

F-^{l-v) ■ 

Iterating, 

p-i fi u \ 

V (1+^)-; ^^™^ foraUrneN, 



F-'^{l-v) 
and so in fact, for any W2 < wi < V, 



F-^l-vi) 
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Putting ui — 1 — vi, U2 = I — V2 and U{C, e) = I — V, we have that for U{C, e) < ui < U2, 

F-\u2) > t\-^°^'^' T^J F~\ui). 
Now whenever z > 1, then [z\ > z/2, so restricting to (1 — ui) > (1 + e)(l — U2) we get 

F-\u2) > 1^'°'"'+''^^ F-\ui), 
which rearranges to the desired result. D 

Lemma 6.4 Fix e > and C > 0. // U2 G (0, 1) is sufficiently close to 1, then for all Ui 
with < ui < (1 + €)u2 — e, 

F-\u,) < F-\u2)il - «l)-^(l - U2f. 

Proof: From Lemma 16.31 we aheady know that this is true when ui > U{C,e). Now if 
ui < U{C, e), then F^^{ui) < F^^{U{C, e)). So it will suffice to show that for U2 sufficiently 
close to 1, the RHS is always at least F~^{U{C, e)). In fact, we will show that the RHS tends 
to 00 as M2 T Ij uniformly in ui. 

For any U2, the RHS is minimised by ui — 0. So we wish to show that 

F-Hu2)il - U2f ^ 00 a.s U2 ^ 1. (6.4) 

Fix wi > C/(2C, e). Then Lemma [6.31 gives . for U2 sufficiently close to 1, 

F-\u2) > F-i(«i)(l - uif^il ~ U2)-2^ 

= c(l-«2)-^^ 

for some constant c. This gives the desired convergence to cxd in (|6.4|l . D 

Lemma 6.5 Let the r.v. U be uniform on (0,1). 

EiF-^{U)\U<u) 

lim , ^ , / — = (J. 

u^i F-i(u)(l-u) 

Proof: Take any C > and e > 0. If m is close enough to 1, then using Lemma [6.41 

E(F-\U);U <u)^ F~\v)dv 
Jo 

p(l-\-e)u — e pu 

< / F-'^{u){l-v)-'^{l-ufdv+ / F-'^(u)dv 

Jo J(l+e)ji-£ 

C+l 

e 



f-^(^)(i-^){ ^'^'r/ + 



We can now choose e as small as desired, and then C as large as desired, to give an upper 

bound on 

E(F-\U);U <u) 

lim sup ., , , , r — 

„^/ F-i(u)(1-m) 
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which is arbitrarily close to 0. 
Finally, for u > 1/2, 

E{F-^{U)\U <u) <2¥.{F-\U);U <u) 

and the result follows. D 

Proof of Proposition l672l We use the representation 

(m|"\...,m(^)) = (f-1([/("\..., [/(:))), 

where ([/{;", . . . , U^^ ) are the order statistics of an i.i.d. sample of size n^ from the uniform 
distribution on (0, 1), written in decreasing order. 
We need to show that 

E F''(U(-'>) > F-i ([/(")) I ^ as n ^ oo. 

Let R > and suppose that u > 1 — R/n? . Then 

n2E(F-i(C/)|C/<u) RK{F-^{U)\U <u) 



F-^u) 



< 



F-^{u)R/n^ 

RE{F-^{U)\U <u) 

F-i(w)(l-u) 

E{F-^{U)\U < u) 



< R sup 

u>l-R/n^ F \U)[1 -U) 



(6.5) 



Now for (almost) all u. 



/n— r+1 



SO for (almost) all il> 1 — R/n? we have, from Markov's inequality, 



E F-\ul:^)>F-\Ml-)) 



L m— r+1 



Jjin) 






t/i") = 



< n- 



,E {F-\U)\U <u) 



< R sup 



i^-i(w) 

E{F-^{U)\U <u) 



>l-R/n^ F-1(u)(1-u) 



using H6.5() . 

So in fact, integrating over u > 1 — R/ri^, 
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E F-\ui:^)>F-\A4"^) 



i'm^r-\-l 



t/,(") >l-R/r^ 



E(F-\U)\U < u) 
< R sup — _. —-^. 



Then 



limsupP( y^ F-i(C/,(,^))>F-i(C/^")) 

'&{F-^{U)\U < u) 



rn— r+l 



< limsupP [t/^") < 1 - R/n^'\ + limsupi? sup 



u>l-R/n^ F H")(l-") 

The second term is by Lemma l^^ for every R. The first term is the probabiUty P(i3„2^^/„2 < 
r — 1), where i3„2 jjy„2 is a binomial {n? ,R/ii?) random variable. This converges as n ^ cxd 
to the probability P(Pr < r — 1) where Pr is a Poisson mean R random variable. Thus this 
probability can be made as small as desired by choosing R large, and hence the limsup is in 
fact as required. D 

6.2 The properties of the greedy path 

In this section we discuss the properties of the greedy path which we have obtained as a 
distributional limit of the optimal path for the discrete problem when a = 0. The path can 
be regarded as a function y = G{x) from [0, 1] to itself in a natural way (specifically, one 
could define G{x) = sup{y ; Yi = {x',y) for some x' < x and some i G A*}. See Figure I^TI 
for a realisation of the greedy path). This function G is monotone non-decreasing, and hence 
defines a measure fi on [0, 1]. We will show that fi is a random self-similar measure which is 
singular with respect to Lebesgue measure, and we will be able to compute its multifractal 
spectrum. 

We recall the definition of the multifractal spectrum for our setting. Let Br (x) denote the 
ball of radius r around the point x S M. For a > 0, we define Ea, the set of points at which 
the measure has local dimension a by 

I rio logr 
The multifractal spectrum is then defined to be 

/(a) = dim_ff (^a), 

where dim/f (A) denotes the Hausdorff dimension of a set A. In our setting we have 

J, / ^rnil r ^ogiGjx + r) - G{x - r)) 

Fa — < X e 0, 1 : lim = a 

[ rio logr 

In particular we note that if fi had a density with respect to Lebesgue measure, then the 
spectrum would be the function /(a) = for all a 7^ 1 and /(I) = 1. 

There are a number of papers making rigorous the multifractal formalism, the heuristic 
argument for computing the multifractal spectrum in terms of the Legendre transform of 
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Figure 6.1: A simulation from the distribution of the "greedy path" which occurs as the hmit 
of the distribution of optimal paths in the case a = 0. 

the moment measures, and we will be able to set our measure in a framework within which 
we can apply this formalism. The study of the multifractal spectrum for random self-similar 
measures is the topic of |lll l^l5] where the underlying assumptions are successively weakened. 
A scaling law on a space E consists of a probability space (fi, !F, P) and for each uj G H. a 
collection of weights and maps {(j)i{uj),pi{uj), . . . ,(/)jv(w),PAr(i^)), where p^ G ]R_|_ and 4>i : E ^ 
E is a contraction with Lipschitz constant r^. For a given scaling law a random self-similar 
measure is a measure fj, which satisfies the distributional equality 



N 



fJ- 



(•) = Ep»^»(^r'(-)), 



i=l 



where Hi are i.i.d. copies of /i (independent of the weights and maps). The support of the 
measure is typically a random self-similar set. 

The multifractal formalism enables the multifractal spectrum for the random self-similar 
measure to be calculated in the following way. Let m{q,9) = E(Ylnp1rf) and j3{q) = \ni{9 : 
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"^(95^) ^ !}• Under the formalism the multifractal spectrum is the Legendre transform of 

f{a)= inUaq + Piq)}. (6.6) 



We will now give a more formal version. 

We introduce a little notation. We write 7^ — {!,..., iV}" for the sequences which 
index the sets after n applications of the scaling law, and T for the tree U^qT^. Let 
(ri®-^,J^® ,P® ) denote the product probability space for random variables on the tree; 
for each node, we have an independent copy of the scaling law. Now for each i S 7^, define 

and 

n = n, (t^0)rj2 (o), J . . . r,„ (wjiJ2...j„_i ). 

The total mass of the random measure over the unit interval is given hyW = liin„^oo X^ier Pi- 
In our setting we will consider random probability measures, so that W = 1. The other limit 
random variable we need is is W{q) = lim„^oo Wn{q) where Wn{q) = X^ier Pf'^i • 

We also define the set Ip C R as follows: q* G Ip if, for some a > 0, the infimum 
infggRJag + /3(g)} is non-negative and is achieved at q*. 

Finally the strong open set condition is that there is an open set O such that, with 
probability 1, one has that 0i(O) n (j)j{0) = for i 7^ j, that ufL^4ii{0) C O and that 

nio) > 0. 

We now state a version of the main result of 8 which can be applied in our setting. 

Lemma 6.6 Let n be a random self-similar probability measure satisfying the strong open set 
condition. If the following three sets of conditions are satisfied: 

1. -c» < EX, .Pi log Pi < 0; 

2. For all q G Ip, EX,(logP,)p?rf ^"^ < 00 and ¥.Y.,{\ogn)plr^^''^ < 00; 

3. For all q £ Ip, E^, ((logp,)^ + (logr,)') pfrf^'^ < 00 and EW^(g) log+ W{q) < 00; 
then the multifractal formalism holds in the following sense. Define 

(3*(a)= inUaq + Piq)}. 



Then for any a > 0, 

/(a)=max{/3*(a),0}, 

with probability 1, and Ea — $ with probability 1 if f3*(a) < 0. 

We now return to the measure arising from our greedy path. From the definition, we can 
describe the greedy path recursively as follows: choose a point Y = (Y{1),Y{2)) uniformly 
in the box [0, 1]^. Then the original path is the union of two independent greedy paths, one 
scaled to lie in [0,^(1)] x [0, r(2)] and the other to lie in [F(l), 1] x [F(2),l]. From this 
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representation we can regard the induced measure as a random self-similar measure for a 
scaling law. Let V and V be independent uniform random variables. Then the scaling law 
has N = 2 with the set of weights {pi,P2) = {V, 1 — V) and the set of contractions {(j)i, (t>2): 
where (j)i{x) — Vx has contraction factor ri = V and 02 (a^) = 1 — (1 — V)x has contraction 
factor r2 = 1 — V . Then the random self-similar measure fi satisfies ^{■) = X]i=i Pil^i{4>^ (')) 
in distribution. From the construction it is clear that the measure is a probability measure 
whose support is the unit interval. 

We note that our measure does not fit into the framework of ^Jj or |2]. The interval 
does not satisfy the strong separation condition as required in ITTl (for a definition see lllp 
but instead it satisfies the strong open set condition with O = (0, 1). This ensures that the 
overlap between a pair of contractions applied to the unit interval occurs at one point. It 
also uses uniform random variables for the contraction ratios and therefore there is no strictly 
positive lower bound on the contraction ratios as required in |2]. 

Theorem 6.7 The measure corresponding to the greedy path has for, any given a > 0, with 

probability 1, 

f\/8^-a-l, 3-2^/2 <a< 3 + 2\/2, 

•^(") =1 n ,h ■ 

I (J, otherwise. 

Ifae [0, 3 - 2\/2) U (3 + 2\/2, cxj) then Ea^% with probability 1. 

Proof: We determine the multifractal spectrum using the multifractal formalism. We need 
to consider m{q,e) = EiVW^ + (1 - ^)'(1 - V)^) and P{q) = M{e : m{q,e) < 1}. It is 
straightforward to compute these quantities and we have 

miq,e) = 2 f y'idy f x'dx= f for g > -1, g > -1. 

Jo Jo (1 + 9)(1 + ^) 



Thus 



and we can calculate that 



m = ^ - 1> 

q+l 



(3* (a) = V8^-a-l 

for all a > 0. We have that /3*(a) is non-negative on the interval [3 — 2-^2,3 + 2\/2] and 
negative elsewhere. Hence, if we can establish the three conditions of Lemma 16.61 to justify 
the formalism, we will have proved our theorem. 

Another calculation gives that If^ — [1 — \/2, 1 -I- \/2] C (—1/2, 3), and that for all q G Ifj, 
one also has f3{q) G (—1/2,3). 

For condition (1), we can compute E^^pilogpi = —1/2. 

For (2) straightforward calculations give E ^j(log pi)p?rf = 2{{q + 1)^(1 + /3))~-^ < oo 
and EX;,(fogn)pfrf = 2{{q + 1)(1 + /3)2)-i < oo, for all q>-l,f3> -1. 

Finally for the conditions (3) we have to do some work. It is easy to calculate the first 
condition 

E^ ((logp.)^ + (iog^.)^).^f - jT^kw) + VTWTW < °°' 

for q > -l,/3 > -1. 
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Thus we only have to verify the final condition. 

We begin by observing that all we need is to prove that for some e > 0, ¥.W{qy~^'^ < oo 
for each q. To do this we will show that EWn{q)^^'^ converges. First we observe that Wn{q) 
is a martingale by the definition of (3{q) and we compute the bracket process 

n 



2 






^i-l 



where Xj ~ V^V + (1 — V,)''(l — Vj)^^'^^ and the (V,, Vj) are independent of each other 
and independent over j. By definition of /3 we know that Exj = 1 and we can also compute 

^■' (1 + 2g)(l + 2/3) ^ r(2g + 2) r(2/3 + 2) 

One can easily check that this quantity is finite over q and (3 in (—1/2, 3), and hence for all 
q (z Ifj. Thus we have 

n 

[W{q)]n =Y,W^-l{2q,2P)(Ex' - 1), 
i=l 

where we write Wn{a, b) = X^ieT Pf'^'t- 

With the bracket process we can control the moments of the martingale as for any 7 > 1 
there is a constant c^ such that 

KWniqy < C,E[W{q)rJ\ 

Thus we need to compute the moments of the bracket process. As we will take 1 < 7 < 2 
and all terms in the sum are positive, straightforward estimates give 

n 

nW{q)]7/' = E(^W^,_i(2g,2/3)(Ex2- 1)^/2 

n 

< E^W,_i(2g,2/3r/2(Ex'-l)T/2 

n 

< ^EW^,_i(7g,7/3)(Ex2- 1)^/2. 
1=1 

Thus we just need to find EW^{-/q,JP) = {EiV'^Wf^ + (1 - V')''^! - V)'^'^))'. Integration 
gives 

Thus our result will hold if we can establish that over the range of q, we can find a 7 > 1 

such that 

2(1 + <Z) ^^ 



{1 + jq)il + q + -f{l - q)) 
A numerical calculation with the quadratic formula shows that this is the case and hence we 
have our result. D 
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The explicit form of the spectrum shows that /(a) > for 3 — 2\/2 < a < 3 + 2\/2- We 
also observe that the set of points for which a = 2 has full dimension 1 . The next result shows 
that the measure corresponding to the path is "singular" : 

Corollary 6.8 With probability 1, the greedy path is continuous and strictly increasing, with 
zero derivative almost everywhere. 

Proof: These properties can be proved fairly directly from the construction of the greedy 
path, but here we deduce them immediately from the multifractal spectrum. 

By construction, the distribution of the path is symmetric in the x and y coordinates; hence 
the continuity property and the property that the path is strictly increasing are equivalent. 

Any point of discontinuity of the path belongs to the set Eq, by definition. But from 
Theorem 16.71 we have that Eq is a.s. empty, so the path is a.s. continuous as desired. 

For the derivative, note that as G is a distribution function it is almost everywhere dif- 
ferentiable with non-negative derivative. Let Di denote the set of points where the path is 
differentiable and has strictly positive derivative. It is straightforward to see that Di C Ei. 
Thus as /(I) = V& — 2 < 1, P — a.s., we have that the Lebesgue measure of Di is with 
probability 1. D 



7 Last-passage random fields and an Airy process 
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Figure 7.1: Figures [7.117. 31 show simulations of the heavy-tailed Airy process Ht for three 
different values of a. Here a = 1. 
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Figure 7.2: The heavy-tailed Airy process Ht in the case a — 1.5. 



In this section we consider the extension of the resuhs to the case of a random field. To 
do this we give a Poisson random measure construction of the continuous limit model. 

Let /i denote a Poisson random measure on R^ with intensity measure \{dxdydz) — 
dxdyaz~°'~^dz. That is, ii N{x,y,z) = fi{{0,x) x (0, y) x (z,oo)) = J^ J^ J^ fj,{dx dy dz) 
denotes the number of points in (0,x) x (0,j/) x (z,oo), then this has a Poisson distribution 
with mean KN{x,y, z) — xyz~°', and the number of points in disjoint sets are independent. 

We can now extend our limiting model to this setting. To relate this to our original model 
we can consider the unit square in R_|_ and order the points of the Poisson random measure 
in decreasing order of their z-coordinates and we recover the sequence of weights Zi — Mi in 
the original model. 

We will write {Yi, Zi) for a point of the Poisson random measure, where the points are 



labelled as above in that we regard Yi as the location in 
continuous last passage percolation model. Let 



and Zi as the weight in our 



C^y = {A C N such that for all Yi ^ Y, for all i,j & A and Yi G [0, x] x [0, y] for all i G A] 
and define 



T{x, y) = sup ^ Z, 



AeC, 



y ieA 



It is clear by Theorem 12. II that this random variable will exist for each fixed x, y. This can be 
extended to show that the random field {T{x, y),x > 0,y > 0} exists and arises as the limit 
of the last passage model. Let 

T^"\x,y)^ max Vx(i;), 



Treni; 



V^TT 
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Figure 7.3: The heavy-tailed Airy process Ht in the case a — 1.99. 



where 11" is the set of directed paths from (1, 1) to ([nx], [ny]), and let 



f(")(.,,) = i^::^i^. 



Theorem 7.1 The field {T{x,y),x > 0, y > 0} exists almost surely and {T^'^'{x,y),x > 
0, y > 0} -^ {T{x, y),x > 0,y > 0} in the sense of convergence of finite dimensional distribu- 
tions. 

Proof: Given a realisation of the PRM, the random variable T{x, y) is well defined for all x, y, 
and over any finite box [0,Cx] x [0,Cy] we have T{x,y) < T{cx,Cy) whenever x < Cx,y < Cy] 
for any fixed Cx, Cy we have P(T{cx, Cy) < oo) ~ 1 and hence by countable additivity 

P(T{x, y) < oo for all x, y) = 1. 

We now need to establish the convergence of finite dimensional distributions. For one 
particular point we already have the one-dimensional convergence result given in Theorem l2.1l 
Exactly the same couplings between the discrete and continuous problems that we used to 
prove Theorem 12 . II can be applied to extend the result to a finite collection of points within 
any finite box. The size of the box is arbitrary and we obtain convergence of all the finite 
dimensional distributions. D 

We now proceed to define a stationary field on the whole of M?. To do this we observe 
that by simple scaling N{Xx, vy, (Xiy)^'" z) = N(x, y, z) in distribution. In particular we have 



T{\x,vy) = {\vY/°'T{x,y) in distribution. 



(7.1) 
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Now put 



&{ui,U2) = exp - 



(mi + U2) 



T(e"\e"=) 



for all ui,U2 G K. Then we have that for any wi, i;2 G M, the collections {8(mi, U2), ui, U2 G M} 
and {&{ui + vi,U2 + V2),ui,U2 G M} have the same distribution. 

The next two results concern the moments and correlations of this stationary field. 



Proposition 7.2 For all (3 G (0, a) and all u G M^ we have Ee{u)i^ = ET(1, l)'' < 00. 

Proof: This is a continuation of the argument given to establish the existence of T{1, 1) in 
the proof of Lemma [3.11 Recalling the setting of that proof, we have T(l, 1) = ^o < Uq, 
where for fc > 0, we defined Uk = J2'iZk+i Li{Mi — M^+i). Here Li is the largest number of 
the first i locations that can be included in an increasing path; there is a constant c such that 
ELi < c\/i and Eif < ci for all i, and the collections {Mi) and (Li) are independent. 

Since r(l,l) < Uq, it will be enough to show that EC/g < cxd. If a finite collection 
of random variables each have a finite /3th moment, then so does their sum. Hence, since 
Uo — Li[Mi — M2) + • • • + Lk{Mk — Mk+i) + Uk, it's enough to show both of the following: 

(i) for any i, E [L,{M, - M,+i)f < 00; 



(ii) for some fc. 



'.uj: <oo. 



For property (i), note that [L^iM^ - M^+i)]'^ < {iMiY , so it's enough to show that E Aff < 
00. But Ml = Wi , where Wi has exponential distribution with mean 1. Thus EMf = 
Jjf" w~^'°'e~'^dw, which is finite for all (3 < a as required. 

So it remains to show (ii). Since /? < a < 2, it is enough to show that EC/^ < 00 for all 
large enough fc, and this is what we will do. 

Recall that we can write Ali ~ {Wi + • • • + Wi)~^^" where Wj are i.i.d. exponential 



Wi, so that Mi 



random variables with mean 1. We also write Vi — W\ + • 

Then, using the fact that (1 + a;)~^/" > 1 — xja for all a; > 0, we have that for all i. 



Y. 



-l/a 



E(Af, 



Mz+i)' 



E 



V, 



-l/a 






-1/a^ 



= E 



<E 



V, 



-2/a 



1-1 



l^H 



V, 



-1/aN 



^rV-vr-Wl,, 



C¥.V, 



-2~2/a 



for some constant C, since Vi and W^+i are independent. Arguing as at (|3.5I) and (|3.6I) . we 
have that Vi has Gamma(i, 1) distribution, and we obtain that for some constant C and all 
large enough i 



E {Mi - M, 



i+i) 



< «-2-2/". 



(7.2) 
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Suppose k is large enough that (|7.2|) holds for all i > k. Then using Cauchy-Schwarz, we 
obtain that for all e > 0, 






i=k+l 



< 



-,1/2 



i=k+l J U=fc+1 

The first sum is finite for any e > 0. Thus squaring and taking expectations, we have that 



i=k+l 



^^"E 



^l+e^.-2-2/a 



oo 
i=*:+l 

for some constants c', c". Since a < 2, this sum is finite for small enough e, and so E C/^ < oo 
as desired. D 



Proposition 7.3 For all (3 < a and all u, v G 



have 



E|r(u)-T(v)|''<4max{||u||,||v||}'3/"||u-v||^/"ET(l)'5. 

Proof: Recall that T(u) is increasing in the partial order on M^. We just need to consider 
the two cases where u < v and where they are not comparable, so that, say, ui < vi and 

U2 > V2- 

For u < V, it is a simple observation that 



T(v)-r(u)<T((0,U2),v)+T((Hi,0),v), 



(7.3) 



where, for x < y G M:^, T(x, y) denotes the maximal weight of an increasing path from x to 
y. By the scaling in the field wc have the distributional relationships 



T{{a, b), (c, d)) ^T{c~a,d~b) = {c- a)i/"(d - 6)i/"r(l, 1). 



(7.4) 



For any random variables Xi, X2, we have E (Xi+Xa)^ < 2'"^'^^''^i> maxjEXf ,EX|}. Thus 
from (|7.3|l and (|7.4|l we obtain 

E|r(v) - r(u)p < 2'"^'^'f'5'i> max {ET(i'i - mi, W2)'', ET(wi, W2 - U2f} 

< 4ET{lf max |(t;i - ui)^/"wf ", wf/"(?;2 - U2)''/"} ■ (7.5) 

For the case ui < vi, V2 < M2, we observe similarly that 



|T(v)-r(u)|<T((ui,0),(t;i,U2))+T((0,«2),(t'i,«2)). 



(7.6) 
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Raising to the power [3 and taking expectations as above, we obtain that 

E |r(v) - r(u)|'' < 4Er(l)^max |(z;i - mf'^v^''', {u2 - i^z)'^/""?^"} • (7.7) 

Combining the estimates from H7.5|l and H7.7|l now gives the resuh. D 

We are now ready to define our analogue of the Airy process. If we set iJ„ = 6('u, —u) = 
T(e",e~"), we have a one-dimensional stationary process {ff„,u G M.}, as the processes 
{Hu+t,u e M) and {Hu,u G M) have the same distribution for all i e M. We note that as 
Hq ~ r(l, 1), the marginal distribution for our stationary process is the distribution of the 
limit random variable in our continuous model for heavy-tailed last passage percolation. By 
applying the estimates for the random field we have estimates on the Holder continuity of the 
heavy-tailed Airy process. 

Corollary 7.4 For (3 < a, we have: 

(i) EH(^ < c». 
(ii) For each r > and all u,v £ [— t, r], 

Proof: The first part is Proposition l7.2l The second follows from the second part of the proof of 
Proposition l7.3l as. assuming u < v and using H7.6|l with vi = e", V2 = e^^, ui = e", U2 — e^", 

E\H,,-K,f < (e''(e-"-e-''))'^/"Ei/^-f ((e''-e")e~")''/"Ei/^ 
= 2(6"-" - l)^/"EiJo^ 
< 2{v - M)^/"e("-")'3/"Eif^, 

and the result follows. D 

Finally we give a weak convergence result for our heavy-tailed Airy process. 

Theorem 7.5 The sequence {a^jT' (e",e~")}tiG[-T,T] converges weakly to {Hu}ue[-T,T] ^'^ 

D[-t,t]. 

Proof: We follow the approach outlined in ,23.^ for such weak convergence problems, in par- 
ticular the proof of ESj Proposition 3.4. 

Since we restrict to u e [— t, r], we can consider only those points of the PRM whose 
locations fall in the box [0, Cj^] x [0,Cy] where c^ = Cy = e'^ . Thus we can write the points 
as a sequence {Yi, Zi),i — 1,2,... such that the sequence (Zi) is decreasing and such that 
Y^ € [0, Cx] X [0, Cy] for all i. 

Define H^ = rfe(e",e-") where 



Tk{x,y) = sup V^i, 

with 

C^y = {^ C {1, . . . , fc} : K, ^ r, Vi, J G A}. 
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As before, we need a corresponding formulation of the discrete model. For given n, we 
work on the box B" = {1, . . . , [ne"^]}^. The weight at a point y E _B" has distribution F, 
independently for different points. We represent the weights and their positions by a vector 
{Yj^,M^, 1 <i < [ne'^]^), where M" form a decreasing sequence. The interpretation is that 
1^" is the location of the ?th largest weight M" in the box B". 

We can then define 77^")'*' = 2\^"^([ne"], [ne~"]) where 



Tt\\'^x^^\'rly^)= sup 



AGCi 






with 



Cif ^{Ac{l,...,kAn^}: Vz, j G A, F/") ~ f/"'}. 



For the proof we will establish that //(")■*-' converges weakly to H'' , that H'' -^ H locally 
uniformly and that 

lim limsupP(p(ff(")'^iJ("') > e) = 0, 

k ^oo ^-i — ^oo 

for each e > 0, where p is the metric for the Skorohod topology. 

We begin by establishing H^ ^- H ask ^ oo. This is a consequence of the construction via 
a PRM. For each u, H^ is an increasing function of k and converges to Hu- With probability 
1, this holds uniformly for each u E [-~t,t], since for all such u we have 



o<H^-m< sup 



E ^^- 



ieA,i>k 



the upper bound is finite with probability 1 exactly as in Lemma l3. II 

Next we wish to show the weak convergence of _ff(")'''' to H''. This can be done by an 
extension of the method of Proposition 13.21 Note that with probability 1 , no two points of 
the PRM share a vertical coordinate or a horizontal coordinate, and in addition no point of 
the PRM falls on the line parametrised by {x, y) = (e", e~"). Then under the same couplings 
used in the proof of Proposition 13. 21 one obtains that with probability 1, ij(")''^' -^ H^ in the 
Skorohod space. Thus this weak convergence also holds as desired. 
Finally we need to control /9(iJ'^"^''^, ij'"^). Fix an e > 0. Consider 



p{H 



(n),k j^(n) 



> e 



< 



sup 

y—T<U<T 



sup 

-T<U<T 



sup 

-T<U<T 



H, 



(n),k 



H, 



(n) 



> e 



f^^"^(e",e"")-f'("'(e",e-") 



sup 



< 



sup 

y—T<U<T 






E^'M^ 



> e 



sup 



' -u «e^ 



Aec' 



T.'^-JM- 



_,, i£A 



> e 



-An) 



where Sj^ '{e^,e ") = sup 



AeC: 



5r(e",e-") 
J2ieA,i>k "n2^*^"- By monotonicity we have 



> e 



sup 

-T<U<T 



stHe^e--) <stHe\en, 
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and using the scaling 

p/p(ij(«),fe^ij(«)^ > e) < p(^si"\l,l) > ee^2r/ay 

By Proposition 13 .31 we have that P(S'^" > ee"^"^/") — > as fc ^ oo uniformly in n. Thus 
indeed we have 

lim limsupP(p(iJ(")^^i/(")) > e) == 0. 

Putting the three pieces together we have shown the weak convergence. D 



8 Higher-dimensional heavy-tailed last passage percola- 
tion 

Up to this point we have considered only two-dimensional models. In this section we indicate 
how to extend most of the results to higher dimensions in a natural way. 

For general d > 2, we consider the passage time from the point (1,1,...,!) to the point 
{n,n,...,n). 

We now consider a sequence of locations i^j , i — 1,2, ... jn"^ which form a uniform 
random permutation of the set {1/n, 2/n, . . . , 1}'^ C [0, 1]'', and a corresponding sequence of 
weights M^', i — 1,2, . . . ,n'^ which are given by the order statistics, in decreasing order, of 
a sample of size n'^ from the underlying weight distribution F . We now assume that the tail 
of F is regularly varying with index a < d. 

Defining C^") as before, we set T^") = sup^gcC") EieA ^^f "^ and f (") = fl-^r^"). 

The continuous model is defined just as before; the locations Yi are now drawn i.i.d. and 
uniformly at random from the box [0, l]** rather than the square [0, 1]^. 

Then T*^") — > T in distribution as n — + cxd; the method of proof is essentially identical to 
that used for Theorem 12. II in the case d = 2. 

The multivariate extensions described in Sections 12.41 and [7| go through in an analogous 
way. For example, we can now obtain a process which is stationary on M"^ such that 



Ul-\ hUrf 

exp ■ 



^-^iy(„) (g«i^ ^ . ,e"^) ,u e R'' I ^ {e(wi, . . .,ua),ue R""} 



as 71 — + cx), in the sense of convergence of finite-dimensional distributions; here T^^' (ui , . . . , Ud) 
is the maximal weight of a path from (1, . . . , 1) to the point ([nwi] , . . . , [nw^]). 

We turn to the path convergence as developed in Section^ ProDOsition l4. II and Theorem 
14. 21 extend easily, with the same method of proof. However, extending Theorem 14. 41 concern- 
ing the convergence of optimal paths viewed as random subsets of [0, l]'^, is more problematic. 
Again we are unable to prove that the optimal path for the continuous model (i.e. the closure 
of UieA* ^) i^ connected (although we expect this to be true). In the case d — 2 this caused 
a little inconvenience but we could work around it by observing that any "jumps" in the path 
consist of horizontal or vertical line segments, and hence that at least there exists a unique 
connected increasing path that contains the optimal path. 
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For d > 3, however, a jump could, for example, cross a square of zero volume in R*^ but 
with non-zero area. Then there is no longer a unique way to extend the optimal path to a 
connected increasing path, and thus there is an ambiguity in the limit object. If we could 
prove the conjecture that the optimal path itself is connected, the convergence in distribution 
of the discrete optimal paths p(")* would follow as before. 

In the case a = 0, we can in fact prove the connectedness of the optimal path for the 
continuous model (defined as in Section using the "greedy algorithm" ) . It is not clear how 
to extend the multifractal analysis of Section f(). 21 However, by analysing a branching random 
walk associated with the algorithm which constructs the greedy path, one can obtain that the 
function from, say, xi € [0, 1] to {x2, ■ ■ ■ , Xd) £ [0, 1]*^^^ which describes the path is almost 
surely everywhere continuous and strictly increasing (although a.s. it also has derivative 
almost everywhere). Thus one can show that p(")* -^ p* in distribution (under the Hausdorff 
metric on subsets of [0, l]'^) for all d in the case a = 0. 
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